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ABSTRACT 

An  extensive  database  on  spatial  structure  in  the  infrared  radiance  of  the  middle  and  upper  atmosphere  has  been 
collected  by  the  Mid-Course  Space  Experiment  (MSX).  The  observed  radiance  contains  spatial  structure  down  to  the 
scale  of  hundreds  of  meters.  This  spatial  structure  results  from  local  fluctuations  in  the  temperature  and  densities  of  the 
radiating  states  of  the  emitting  molecular  species  as  well  as  fluctuations  in  radiation  transport  from  the  emitting  regions 
to  the  observer.  A  portion  of  this  database  has  been  analyzed  to  obtain  statistical  parameters  characterizing  stochastic 
spatial  structure  in  the  observed  radiance.  Using  simple  models,  the  observed  statistics  have  been  shown  to  agree  with 
prior  observations  and  theoretical  models  of  stochastic  spatial  structure  generated  by  gravity  waves  for  special  viewing 
geometries.  The  SHARC  model  has  been  extended  to  predict  the  statistics  of  stochastic  fluctuations  in  infrared  radiance 
from  the  statistics  characterizing  temperature  fluctuations  in  the  middle  and  upper  atmosphere  for  arbitrary  viewing 
geometries.  SHARC  model  predictions  have  been  compared  with  MSX  data  and  shown  to  be  in  generally  good 
agreement.  Additional  work  is  in  progress  to  account  for  the  statistics  characterizing  small  spatial  scale  fluctuations. 

Keywords:  Middle  and  upper  atmosphere,  infrared,  spatial  structure,  gravity  waves,  Mid-Course  Space  Experiment, 
SHARC 


1.  INTRODUCTION 

Gravity  waves  play  a  critical  role  in  determining  the  structure  of  the  lower  and  middle  atmosphere  and  also  play  a  key 
role  in  determining  the  spatial  structure  of  the  radiance  of  the  upper  atmosphere.  Small  scale  waves  are  a  dominant 
source  of  the  transport  of  atmospheric  energy  and  momentum  A  new  and  important  data  base  on  gravity  waves  in  the 
middle  and  upper  atmosphere  has  been  collected  by  the  Mid-Course  Space  Experiment  (MSX)3,4.  The  first  space-based 
mid-wave  infrared  (MWIR)  observations  of  gravity  waves  were  made  using  the  narrow  band  MWIR  sensor  on  MSX5. 
Circular  gravity  waves  generated  by  a  thunderstorm  cell  have  also  been  observed  by  MSX6.  The  extensive  MSX  data 
base  on  spatial  structure  in  atmospheric  radiance  has  now  been  analyzed  to  determine  statistics  characterizing  stochastic 
spatial  structure  in  the  radiance  of  the  middle  and  upper  atmosphere.  Theoretical  approaches  for  relating  the  power 
spectral  densities  (PSDs)  of  the  observed  radiance  fluctuations  to  the  underlying  temperature  and  density  fluctuations 
have  been  developed.  The  inferred  temperature  statistics  have  been  compared  both  to  prior  measurements  and 
theoretical  accounts  of  the  statistics  of  temperature  and  wind  velocity  fluctuations  caused  by  gravity  waves  in  the 
stratosphere  and  mesosphere. 

2.  MSX  OBSERVATIONS  OF  SPATIAL  STRUCTURE  IN  MIDDLE  AND  UPPER 

ATMOSPHERIC  RADIANCE 

2.1  The  Mid-Course  Space  Experiment  (MSX) 

The  MSX  satellite  was  launched  into  a  900  km  altitude,  nearly  sun-synchronous,  circular  polar  orbit.  A  key  mission 
goal  was  collecting  data  on  the  spatial  structure  in  the  radiance  of  the  middle  and  upper  atmosphere  using  the  SPIRIT-3 
radiometer.  The  six  wavelength  channels  of  the  SPIRIT-3  sensor  included  the  B1  band  (4.22-4.36  pm)  which  was 
selected  to  collect  data  on  atmospheric  radiance  in  the  heart  of  the  4.3  pm  C02  absorption  band.  The  contributions  to 
radiance  observed  in  the  B1  band  come  from  altitudes  above  20  km  due  to  optical  opacity  effects.  B1  band  data  was 
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collected  for  a  wide  variety  of  viewing  geometries  including  below  the  horizon  (BTH)  near-nadir  and  near  limb  as  well 
as  above  the  horizon  (ATH)  viewing  geometries.  Sensor  saturation  effects  limited  data  collected  on  atmospheric 
radiance  from  03  in  the  SPIRIT-3  A  band  (6.8-10.8  Jim)  to  ATH  viewing  geometries  with  tangent  heights  above  55  km. 
Each  radiometer  channel  consisted  of  two  or  more  staggered  linear  arrays  of  192  pixels  each  having  a  field  of  view  of  90 
prad  (the  B  band  was  divided  into  B1  and  B2  bands  each  consisting  of  two  rows  of  approximately  86  pixels, 
respectively).  The  spatial  resolution  was  80  and  300  m  per  pixel  for  nadir  and  earth  limb  views,  respectively.  The  data 
presented  here  was  collected  in  a  push  broom  scan  mode  with  the  detector  arrays  oriented  perpendicular  to  the  spacecraft 
velocity  as  illustrated  in  Fig.l.  Great  care  was  put  into  pre-launch  and  on-orbit  calibration  of  the  SPIRIT-3  sensors,  and 
the  uncertainty  in  the  absolute  radiances  measured  by  the  SPIRIT-3  radiometer  is  less  than  10%.  IR  data  collections 
were  conducted  using  the  cryogenic  SPIRIT-3  sensors  on  MSX  from  24  April  1996  until  25  February  1997. 

2.2  MSX  data  collections  used 

Data  presented  in  this  paper  was  collected  using  the  MSX  SPIRIT-3  B1  and  A  bands  for  a  variety  of  viewing  geometries. 
BTH  and  ATH  data  collected  in  the  narrow  MWIR  B1  band  was  analyzed  to  determine  the  effects  of  temperature  and 
density  fluctuations  on  radiance  from  C02.  ATH  data  collected  in  the  long  wave  infrared  (LWIR)  A  band  was  also 
analyzed  to  determine  temperature  and  density  effects  on  03  radiance  fluctuations  for  comparison  purposes.  For  this 
analysis  the  following  BTH  scene  data  was  used: 


Band 

Nadir  Angle 

Scene  Dimensions 

Scene  1 

B1 

1.3° 

6.7  x  2701  km2 

Scene  2 

B1 

55° 

13.5  x  1164  km2 

Scene  3 

B1 

57° 

15.1  x  4237  km2 

In  addition,  the  following  ATH  scene  data  was  analyzed: 


Band 

Tangent  Height 

Scene  Dimensions 

Scene  4 

B1 

1.8-16  km 

14.2  x  1050  km2 

Scene  5 

B1 

40  -  54  km 

14.0  x  1367  km2 

Scene  5A 

A 

57.3 -71.2  km 

13.9  x  1063  km2 

The  MSX  line  of  sight  and  tangent  point  ground  tracks  for  these  data  collections  are  shown  in  Fig  2.  Shown  in  Fig.  3  are 
the  six  MSX  scenes  used  for  the  present  analysis  of  atmospheric  structure. 

2.3  MSX  scene  data  analysis 

The  six  MSX  scenes  described  above  were  analyzed  to  obtain  a  variety  of  descriptors  of  the  spatial  content  of  the  scenes. 
The  scene  data  was  sufficiently  over-sampled  during  collection  to  maintain  the  spatial  resolution  inherent  in  the  sensor 
footprint  and  to  avoid  aliasing  effects.  Average  radiance  and  relative  standard  deviation  profiles  derived  from  the  pixel 
time  histories  are  shown  in  Fig.  4.  The  relative  standard  deviation  profiles  include  both  deterministic  and  gravity  wave 
induced  stochastic  structure  and  thus  set  an  upper  limit  on  the  stochastic  structure  in  the  scenes. 

Data  over-sampling  permitted  decimation  and  co-adding  of  data  samples  in  the  in-scan  direction  to  improve  signal  to 
noise  and  thus  the  power  spectral  densities  (PSDs)  derived  from  the  data.  To  determine  the  statistical  parameters 
characterizing  the  scenes,  each  scene  was  de-trended  by  subtracting  a  first  or  second  order  polynomial  fit  to  the  radiance 
for  each  row  in  the  scene.  Shown  in  Fig.  5  is  a  segment  of  the  de-trended  scenes  presented  with  unity  aspect  ratio  (equal 
distance  scales  in  the  in-scan  and  cross-scan  directions).  As  the  viewing  geometry  changes  from  nadir  viewing  to  limb 
viewing  the  spatial  structure  in  the  de-trended  scenes  changes  from  isotropic  to  striated. 

Both  classical  and  modem  spectral  analysis  techniques  were  used  to  estimate  PSD  functions  as  a  function  of  wave 
number  (km1)  for  the  observed  stochastic  radiance  structure.  Averaged  periodograms  of  the  structure  data  were 
generated  using  the  standard  techniques  of  pre-whitening  the  data7,  applying  a  Hanning  data  window8,  and  post-coloring 
the  periodogram9.  Spectral  noise  was  reduced  by  computing  separate  periodograms  for  each  column  and  row  of  each  de- 
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trended  scene  and  averaging  the  sets  of  row  and  column  periodograms,  respectively.  Modem  spectral  analysis 
Dc^™qUes  were  also  used  in  whlch  Burg’s  method8  9  was  used  to  estimate  the  coefficients  of  an  auto-regressive  (AR) 
PSD  function.  Prior  to  computing  PSDs,  a  number  of  steps  were  taken  to  reduce  the  non-stationarity  of  the  data.  Each 
scene  was  de-trended,  as  described  above.  Then  the  de-trended  values  for  each  row  were  divided  by  the  standard 
deviation  for  the  row  and  multiplied  by  the  average  standard  deviation  for  the  entire  scene.  Since  the  data  are  non- 
stationary  in  the  cross-scan  direction,  cross  scan  data  were  limited  to  45  pixels  for  ATH  data  and  41  or  84  pixels  for 
BTH  data.  The  in-scan  and  cross-scan  periodogram  and  auto-regressive  PSDs  for  the  six  scenes  are  shown  in  Figs.  6 
and  7.  The  in-scan  (horizontal)  large  spatial  scale  PSD  slopes  are  ~  -5/3  for  nadir  viewing  geometries  and  -  -  8/3  for 
imb  viewing  geometries.  The  in-scan  PSD  curves  exhibit  a  change  in  slope  at  a  spatial  scale  of  10  km.  The  cross-scan 
PSD  slopes  are  ~  -4  for  limb  viewing  geometries. 

2.4  Consistency  of  MSX  radiance  observations  with  observations  and  prior  theory  of  gravity  waves 

Prior  observations  of  gravity  waves  show  temperature  and  wind  velocity  PSDs  generally  varying  as  k/5/3  10  although 
there  is  evidence  of  kx‘  behavior  at  small  spatial  scales  between  1  and  10  km  under  some  conditions'1-12-13’'4-15'1’6'17'18'19-® 
and  k,  •  ,  where  kx  and  kz  are  wavenumbers  in  horizontal  and  vertical  directions.  A  number  of  theories  of  gravity 

waves  also  predict  temperature  and  wind  velocity  PSDs  varying  as  kx'5/3  and  k/3  However 

the  SPIRIT-3  sensor  on  MSX  measures  radiance  observed  along  lines  of  sight,  not  in-situ  temperature  or  wind  velocity. 
Calculations  of  radiative  contribution  functions  for  atmospheric  4.3  pm  radiance  show  that  the  observed  MWIR  radiance 
m  the  B1  band  comes  from  a  narrow  range  of  altitudes  around  40  km.  The  nadir  to  limb  viewing  observations  made  by 
MSX  span  a  wide  range  of  ratios  of  path  length  in  the  radiating  layer  to  the  spatial  scale  of  atmospheric  temperature 
fluctuations.  A  model  calculation  of  the  effects  of  integrating  through  an  isotropic  two  dimensional  fluctuation  field 
(displayed  in  Fig  8)  shows  that  the  PSD  changes  by  a  factor  of  k'1  when  the  path  length  in  the  radiating  layer  is  on  the 
order  of  2  to  3  correlation  lengths.  The  steepening  of  the  PSD  slope  can  be  understood  as  the  result  of  the  low-pass 
filtering  resulting  from  the  integration  along  the  line  of  sight.  This  result  suggests  that  the  change  in  the  slopes  of  the 
large  spatial  radiance  PSDs  observed  by  MSX  from  -5/3  for  nadir  viewing  geometries  to  -8/3  for  limb  viewing 
geometries  is  a  filtering  effect  due  to  observing  a  fluctuating  layer  of  fixed  thickness  along  a  slant  path.  A  parallel 
argument  applies  to  why  the  slopes  of  the  small  spatial  scale,  cross-scan  (near  vertical  in  the  BTH  near-limb  and  ATH 
scans)  radiance  PSDs  are  observed  to  be  -4,  as  opposed  to  -3  which  would  be  expected  for  the  kz  dependence  of  the 
temperature  spectrum.  However  a  detailed  atmospheric  model  is  required  to  relate  temperature  and  density  fluctuation 
spectra  in  the  middle  and  upper  atmosphere  to  observed  radiance  fluctuation  spectra  for  arbitrary  viewing  conditions. 

3.  SHARC  MODEL  FOR  STOCHASTIC  SPATIAL  STRUCTURE 


3.1  Introduction 

The  Air  Force  Research  Laboratory  originally  developed  the  SHARC  model36  to  describe  line  of  sight  spectra,  in-band 
radiance,  and  transmittance  in  the  middle  and  upper  atmosphere.  A  key  feature  of  the  SHARC  model  is  that  it  describes 
non-local  thermodynamic  equilibrium  (NLTE)  effects  on  the  populations  of  radiating  molecular  vibrational  states,  which 
is  important  for  all  species  above  some  altitude  in  the  middle  and  upper  atmosphere  but  not  important  in  the  lower 
atmosphere.  The  SHARC  model  has  subsequently  been  extended  to  include  stochastic  spatial  structure  due  to 
superposed  and  interacting  atmospheric  gravity  waves37-38.  A  key  challenge  is  developing  a  computationally  tractable 
method  for  relating  temperature  and  density  fluctuations  in  a  sensor  field  of  view  to  the  fluctuations  in  the  radiance  in 
e  image  plane  of  the  sensor.  SHARC  provides  the  capability  for  predicting  scene  spatial  structure  in  atmospheric 
radiance  in  the  wavelength  range  2-40  pm  for  arbitrary  viewing  geometries  as  long  as  there  is  no  significant  contribution 
to  the  radiance  from  altitudes  below  30  km.  Thus  the  SHARC  model  is  appropriate  for  predicting  stochastic  spatial 
staicture  observed  by  MSX  in  the  B1  band  for  the  full  range  of  viewing  geometries  and  for  structure  observed  in  the 
MSX  A  band  for  tangent  heights  above  30  km.  As  will  be  shown  below,  SHARC  model  predictions  of  radiance  PSDs 
are  consistent  with  analytic  solutions  from  gravity  wave  theory  for  special  limiting  cases  (nadir  and  horizontal  viewing 
geometries).  The  model  predictions  are  also  consistent  with  the  MSX  MWIR  observations  of  large-scale  BTH  stochastic 
spatial  structure  and  of  small-scale  stochastic  structure  observed  for  ATH  lines  of  sight. 


The  inputs  to  the  new  SHARC  model  for  stochastic  spatial  structure  are  the  data  base  of  the  altitude  profiles  of  kinetic 
temperature  and  atomic  and  molecular  species  densities  provided  by  the  SHARC  Atmospheric  Generator  (SAG)43  the 
data  base  on  molecular  spectral  line  parameters  provided  by  H1TRAN-9244,  and  a  model  of  three  dimensional  structure 
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describing  the  atmospheric  temperature  fluctuation  spectra  and  autocovariance  functions.  The  PSDs  are  assumed  to  be 
separable  into  horizontal  and  vertical  components,  and  the  altitude  dependence  of  temperature  fluctuation  correlation 
lengths  and  variance  is  provided  by  Strugala  et.  al.39.  The  SHARC  spectral  model  is  used  to  calculate  the  populations  of 
the  vibrational  levels  of  radiating  species  under  NLTE  conditions  and  the  resulting  radiance  spectra  and  in-band  radiance 
as  a  function  of  the  line  of  sight.  The  SHARC  spatial  structure  model  generates  radiance  autocovariance  contribution 
functions,  which  enable  one  to  calculate  radiance  PSDs,  variances  and  correlation  lengths  for  non-stationary  stochastic 
radiance  fluctuations  due  to  gravity  waves. 

3.2  Details  of  the  SHARC  structure  model  calculation 

If  temperature  fluctuations,  AT(r),  are  sufficiently  small  that  linear  response  can  be  assumed,  the 
radiance  fluctuations,  L^p)  ,  in  the  pass  band,  AA.,  in  an  image  plane  can  be  expressed  as  an  integral  along  the  line 
of  sight  of  the  line  of  sight  radiance  response  function,  F^(r) ,  which  accounts  for  all  path -integration  and  geometrical 
filtering  effects,  times  the  atmospheric  kinetic  temperature  fluctuations,  A T(r) : 

(1)  MhX(p)  =  \drF^{r^nr). 

Here,  the  vector,  p ,  defines  the  pixel  location  in  the  image  plane,  and  r  is  the  position  vector  of  the  pixel  p  relative  to 
the  source  point.  It  has  been  assumed  that  dynamical  fluctuations  are  adiabatic,  so  that  density  fluctuations  can  be 
expressed  simply  in  terms  of  temperature  fluctuations.  In  the  case  of  species  and  conditions  where  NLTE  effects  are 
important,  FA^(r)  includes  contributions  from  all  kinetic  temperature  fluctuations,  including  those  which  enter  through 
the  dependence  of  the  vibrational  temperatures  characterizing  the  populations  of  the  radiating  state  on  kinetic 
temperature.  In  general,  it  is  possible  for  the  radiative  transport  of  temperature  fluctuations  in  regions  not  on  the  line  of 
sight  to  produce  additional  radiance  fluctuations  along  the  line  of  sight.  Then  equation  (1)  becomes  a  volume  integral 
rather  than  a  line  integral.  We  will  specifically  exclude  such  nonlocal  effects  in  this  paper,  although  it  is  possible  that 
they  might  play  some  role  in  Scenes  5  and  5A. 

The  radiance  autocovariance  function,  CovL(p,p') ,  can  be  expressed  as  CovL(p,p')  =  £[ALA^(^)ALA^(p')]  where  p 
and  p  point  to  two  locations  in  the  sensor  image  and  E[x]  is  the  expectation  value  of  x.  From  equation  (1), 
CovL(p,p')  can  then  be  expressed  in  terms  of  line  of  sight  integrations  of  the  atmospheric  temperature  autocovariance 
function,  CovT  (r ,  r *) ; 

(2)  CovL  (p,  p')=  j  dr  [J  dr’Fvi  (F)^U  (r)CovT  (r.r')]. 

Here  r  and  7  are  two  line  of  sight  vectors  that  point  to  the  p  and  p  end  points  in  the  image  plane.  The  quantity  in 
square  brackets  is  the  line  of  sight  autocovariance  contribution  function,  WL(r,p, p') ,  so  that 

Covi(p,p')  =  jdrWi(r9p9p') .  Under  conditions  of  local  stationarity  where  the  characteristic  spatial  scales  of  the 

fluctuating  quantities  is  small  compared  to  the  scale  over  which  the  mean  values  of  these  quantities  change  (for  example, 
the  vertical  temperature  correlation  length  is  much  less  than  the  atmospheric  scale  height),  the  expressions  for  the 
radiance  and  temperature  autocovariance  functions  depend  strongly  only  upon  the  displacements  A/3  =  /?-/?'  and 
A r  —  r  —  7. 

Under  conditions  of  local  stationarity,  taking  the  Fourier  transform  of  equation  (2)  leads  to  a  simple  relation  between  the 
radiance  PSD  and  the  three  dimensional  temperature  PSD: 

(3)  PSDL(p,kv,kh)  =  J dr\±j dklose,k‘°*r  FaX (F)QaX (klos )PSDT (F, kloskv,kh )  , 

where  Q^{kios)  is  the  Fourier  transform  of  the  LOS  radiance  response  function,  F^(r'),  kios  is  the  spatial  wave 
number  along  the  LOS,  kv  and  khare  the  “sensor  vertical”  and  “sensor  horizontal”  wave  numbers  respectively,  and  where 
the  sensor  horizontal  direction  is  the  direction  in  the  image  plane  parallel  to  the  horizon  and  the  sensor  vertical  direction 
is  perpendicular  to  both  the  sensor  horizon  and  the  line  of  sight.  If  the  atmospheric  temperature  PSD  is  stationary  over 
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that  portion  of  the  line  of  sight  which  contributes  to  the  radiance  then  the  dependence  of  the  radiance  PSD  on  range 
along  the  line  of  sight  can  be  neglected,  and  equation  (3)  can  be  re-written,  after  integrating  over  r, 


(4) 


PSD L  (*W .  h  )  =  J-  J  dkhs  |  Q&X  (klos )  | 2  PS  Dr  (klos  ,kv,kh). 


Here  | Qax(^Ios)  »  path-integration  filter  in  spatial-frequency  space, 
transform  of  the  line  of  sight  radiance  response  function. 


is  the  square  of  the  magnitude  of  the  Fourier 


As  described  above,  a  key  input  to  the  SFLARC  stochastic  radiance  model  is  the  three  dimensional  temperature 
autocovariance  function.  Assuming  that  the  temperature  autocovariance  function  is  separable  into  horizontal  (h=(x,y)) 
and  vertical  (z)  components,  then  CovT(r,  Ar)  =  CovH(r,  AH)Covz(r,Az) ,  and  the  corresponding  three-dimensional 

PSD  is  also  separable  and  can  be  expressed  as:  PSDT(r,k)  =  PSDH(r,kH)PSD2(r,kz) .  The  von  Karman  spectral 

model40,41,42  is  used  here  to  model  the  atmospheric  vertical  and  horizontal  power  spectra  and  autocovariance  functions. 
The  vertical  temperature  PSD  and  autocovariance  functions  are  given  by: 

(5)  PSDz(r,kz)  =  Nzal  (r)tz  (r  XI + (?)*?  )"(v' +U2) 

and 

(6)  Covz(r,Az)  =  MzO$(r)(Az/tz(r))v’  KVf  (A zUz{f)) . 

Similarly,  the  horizontal  temperature  PSD  and  autocovariance  functions  are  given  by: 

(7)  PSDH(r,kH)=NHa}(r)e]f(n{l  +  e]J(r)k2H)-{v"+  »> 

and 

(8)  CoVff(r ,AH)  =  MHo}(r)(AH UH{r))v "  KVft  (A H/t„(r))  ■ 

where  tyj(r)  =  i2x(r)  +  (r) , kff  =  kx  +  k2  ,  and  AH2  =  Ax2  +  Ay2.  Here  Kv*  and  KVf)  are  modified  Bessel 

functions  of  the  second  kind  of  order  vz  and  vH ,  respectively,  and  NZi  Mzy  NHi  and  MH  are  normalization  constants.  The 
temperature  variance,  <T^(r),  vertical  scale  length,  fz(r),  and  horizontal  scale  length,  depend  parametrically 

on  r  as  the  LOS  continuously  changes  altitude.  Values  for  the  scale  lengths  and  temperature  variances  have  been  taken 
from  Strugala  et  a/39..  By  selecting  the  Bessel  function  orders  vz  and  vH  as  1  and  1/3  respectively,  the  asymptotic  vertical 
spectral  index  at  large  k^  is  -3  and  the  2-D  asymptotic  horizontal  spectral  index  is  -8/3,  leading  to  a  1-D  horizontal 
spectral  index  of -5/3.  The  model  is  consistent  with  simple  gravity  wave  model  predictions  and  observations  for  spatial 
wave  number  regions  where  the  vertical  index  is  -3  and  the  1-D  horizontal  index  is  -5/3.  The  assumption  of  separability 
in  kx  and  k^  is  a  reasonable  approximation  for  the  case  of  multiple  superposed  gravity  waves  as  opposed  to  a  single 
gravity  wave. 


Because  of  the  non-stationarity  and  the  lack  of  vertical-horizontal  isotropy,  observational  effects  on  radiance  PSD  slopes 
due  to  line-of-sight  integration  are  a  function  of  the  relative  vertical  and  horizontal  scale  lengths  and  the  orientation  of 
LOS  through  the  atmospheric  structure.  There  are  two  special  orientations  for  which  analytical  results  can  be  obtained. 
The  first  situation  where  one  can  obtain  an  analytic  solution  is  pure  nadir  or  zenith  viewing,  in  which  the  LOS 
integration  is  in  the  vertical  and  hence  independent  of  horizontal  structure.  There  is  no  steepening  of  the  radiance  PSD 
spectral  slope,  and  nadir  views  through  the  separable  gravity  wave  model  have  a  radiance  PSD  with  a  spectral  slope 
close  to  -5/3.  The  second  situation  where  analytic  solutions  are  possible  is  pure  horizontal  viewing,  in  which  the  LOS 
passes  through  the  isotropic  horizontal  structure.  The  1-D  horizontal  wave  number  spectral  slopes  steepen  with  the  index 
changing  from  -5/3  to  -8/3  38.  Limb  viewing  often  approximates  pure  horizontal  viewing  since  the  limb  path  is 
essentially  a  horizontal  path  at  the  tangent  point.  It  was  shown38  that  for  power  law  type  PSDs,  the  magnitude  of  the 
spectral  index  increases  by  one.  Hence,  for  limb  viewing,  when  the  radiance  band  is  sufficiently  optically  thin  that  the 
contributions  come  primarily  from  the  tangent  point  altitude,  the  vertical  power  spectral  index  is  essentially  unaffected 
and  the  one  dimensional  horizontal  index  changes  from  -5/3  to  -8/3.  The  LOS  integration  also  transforms  the  order  of 
the  modified  Bessel  functions  occurring  in  the  radiance  autocovariance  function.  The  order  is  increased  by  Vi 
corresponding  to  the  spectral  index  changing  from  -n  to  -(/?+l). 
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4.  DATA-MODEL  COMPARISONS 


4.1  Comparisons  of  Radiance  and  Their  Relative  Standard  Deviations 

The  full  SHARC  model  was  used  to  predict  radiances  and  standard  deviations  of  stochastic  fluctuations  in  radiance  for 
the  six  MSX  scenes  described  in  Section  2.2.  Measured  and  predicted  radiances,  as  well  as  relative  measured  and 
predicted  standard  deviations  of  the  stochastic  component  of  radiance  are  shown  in  Table  1.  The  measured  and 
predicted  radiances  agree  within  the  uncertainties  of  the  model  altitude  profiles  of  temperature.  All  discrepancies  are 
accounted  for  by  temperature  adjustments  of  about  5  K  at  40  km.  However,  the  predicted  relative  standard  deviations 
for  radiance  fluctuations  are  significantly  larger  than  the  relative  standard  deviations  obtained  from  the  six  MSX  scenes. 
The  predicted  standard  deviations,  of  course,  reflect  the  temperature  variance  database  in  the  model39.  Analysis  of  the 
extensive  MSX  scene  database  will  provide  additional  inputs  for  adjusting  the  temperature  variance  database. 

4.2  Radiance  PSD  Comparisons 

The  principal  contributions  to  radiance  fluctuations  in  the  narrow  4.3  pm  B1  band  for  the  viewing  geometries  of  Scenes 
1  -  4  come  from  an  atmospheric  layer  with  maximum  contribution  at  40  km.  In  addition,  day  lit  portions  of  viewing 
lines  of  sight  also  have  NLTE  contributions  from  the  65  to  80  km  altitude  range  arising  from  solar-excited  fluorescent 
emissions.  Shown  in  Fig.  9  is  the  LTE  portion  of  the  radiance  variance  contribution  function  calculated  for  the 
conditions  of  Scenes  1  and  3.  The  contribution  function  is  slightly  broader  for  nadir  than  for  sub-limb  viewing. 

When  contributions  from  broad  or  multiple  layers  are  present,  the  non-stationarity  of  the  temperature  and  density 
structure  of  the  atmosphere  needs  to  be  taken  into  account.  The  effect  of  path  integration  is  to  act  as  a  low  pass  filter. 
Thus  the  radiance  PSDs  can  be  expected  to  have  a  steeper  fall  off  with  wave  number  than  do  the  underlying  temperature 
PSDs  for  viewing  geometries  where  the  line  of  sight  passes  obliquely  through  a  radiating  layer  or  lies  nearly  horizontal 
to  it.  The  dependence  of  the  radiance  PSD  spectral  index  on  the  local  zenith  angle  (LZA),  or  angle  that  the  line  of  sight 
makes  with  the  radiating  layer  at  40  km  for  the  in-scan  and  cross-scan  directions,  is  shown  in  Figs.  10  and  1 1  for  two 
ratios  of  horizontal  to  vertical  temperature  correlation  lengths.  The  slope  of  the  PSD  varies  from  -5/3  to  -8/3  in  the  in¬ 
scan  direction  and  from  -5/3  to  -4  in  the  cross-scan  direction,  as  has  been  described  above.  Also  shown  are  the  in-scan 
radiance  spectral  slopes  obtained  from  MSX  scenes  1  to  5 A  for  spatial  scales  >  ~10  km,  which  agree  reasonably  well 
with  the  model  predictions.  The  model  predictions  for  the  slopes  of  the  cross-scan  radiance  PSDs  also  agree  well  with 
the  measured  PSD  slopes  shown  in  Fig.  7. 


5.  CONCLUSIONS 

Unique  data  on  the  spatial  structure  of  the  radiance  of  the  middle  and  upper  atmosphere  has  been  collected  using  the 
SPIRIT-3  radiometer  on  MSX.  The  radiance  fluctuation  observations  made  by  MSX  are  consistent  with  stochastic 
spatial  structure  resulting  from  atmospheric  gravity  waves.  The  SHARC  model  used  to  predict  infrared  spectra,  in-band 
radiance,  and  transmittance  in  the  middle  and  upper  atmosphere  has  been  extended  to  predict  stochastic  spatial  structure 
due  to  gravity  waves.  The  SHARC  predictions  of  stochastic  spatial  structure  agree  relatively  well  with  the  MSX  data  for 
large-scale  spatial  structure  (scale  lengths  >  ~10  km).  A  new  capability  is  the  ability  to  relate  the  variances  and  PSDs  of 
infrared  radiance  in  the  middle  and  upper  atmosphere  for  arbitrary  viewing  geometries  to  the  variances  and  PSDs 
describing  atmospheric  temperature  fluctuations.  With  this  capability  the  extensive  MSX  database  on  spatial  structure  in 
the  atmospheric  radiance  can  be  fully  exploited  to  characterize  the  temperature  structure  of  the  middle  and  upper 
atmosphere.  Additional  work  is  in  progress  to  account  more  fully  for  the  small-scale  spatial  structure  observed  by 
MSX. 
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Scene 

# 

Approx. 

latitude 

start/end 

Local  zenith 
angle 

@  40  km  alt 
(deg.) 

Tan.  Alt. 
(km) 

Measured1 

Radiance 

(Trend) 

(>iW/cm2-sr) 

Predicted  avg. 

radiance 

(fiW/cm2-sr) 

Measured2 
Std.  Dev./ 
Row  Avg. 

% 

Measured3 
Std.  D c\J 
Trend 

% 

Predicted 
Std.  Dev.  / 
Predicted 
Avg. 

% 

1 

1 2°N  -  35°N 

1.7 

BTH/ 

nadir 

.90-  1.06 
.90-1.06 

1.0 -1.2 

1.04 

1.21 

1.11-0.94 

1.29-1.09 

4.9 

2 

30°N  -  40“N 

78 

BTH/ 

terminator 

.77  -  .53 
.74 -.52 

0.94  -  0.78 

1.32 

1.50 

1.11-1.61 

1.28-1.81 

5.1  -5.3 

3 

I9°N  -  55°N 

82 

BTH 

1.12-1.28 
1.13-  L29 

1.3- 1.6 

1.33 

1.40 

1.41-1.23 

1.50-1.32 

4.2 

4 

60°N  -  52°N 

84-85 

2 

16 

1.32-1.20 

1.35-124 

1.7-  1.6 

1.8-  1.7 

0.95 

0.89 

0.92-1.01 

0.86-0.94 

4.2  -  4.4 

5 

15°N  -  3°N 

90 

40 

54 

1.33-1.44 
0.36  -  0.23 

1.8-  1.9 
0.31-0.34 

2.97 

5.32 

3.04-2.81 
11.1  -6.70 

3.6 

8.2 

5A 

1 5°N  -  5°N 

- 

57 

71 

8.67-8.77 

0.94-0.85 

10.2-9.9 
1.5- 1.4 

2.36 

1.33 

2.44-2.44 

1.25-1.39 

4.0 

3.0 

Table  1.  Comparison  of  predicted  versus  measured  average  radiance  and  the  percent  relative  standard  deviation  of  the  radiance 
for  six  MSX  scenes.  Standard  deviations  are  calculated  from  detrended  images. 


1 .  Estimate  of  the  average  radiance  from  trend  at  4  comers  of  the  scene. 

2.  Estimate  of  the  standard  deviation  of  the  radiance  (removing  the  trend)  /  row  average  (%). 

3.  Estimate  of  the  standard  deviation  of  the  radiance  (removing  the  trend)  /  trend  (%)  at  4  comers  of  the  scene. 


Figure  1.  The  projection  of  the  MSX  focal  plane  for  the  MWIR 
bands,  B|  and  B2,  on  the  lower  atmosphere  is  shown  for  a  case 
where  the  line  of  sight  is  in  a  plane  orthogonal  to  the  orbit  and 
elevated  55  degrees  from  the  nadir. 


LONGITUOE 


Figure  2.  The  MSX  line  of  sight  target  region  ground 
track  or  tangent  point  tracks  are  shown  for  the  below 
the  horizon  (BTH)  and  the  above  the  horizon  (ATH) 
events  presented  here.  The  scene  number,  the  local 
zenith  angle  (LZA)  at  the  ground  for  Scenes  1-3  and 
the  tangent  altitude  (tan  ht)  for  Scenes  4-5  and  the  solar 
zenith  angles  (SZA)  are  indicated. 
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Scene  1,  May  25.  1996,  MWIR,  Nadir  1.3  Deg 
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Scene  2.  Nov  13.  1996,  MWIR,  Nadir  55.0  Deg 
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Scene  3,  Aug  17,  1996,  MWIR,  Nadir  57.0  Deg 
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Scene  4,  Aug  19,  1996,  MWIR,  Nadir  61.0  Deg,  Tan  Alt  1.8  -  16.0  (Km) 
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Scene  5,  Aug  10,  1996,  MWIR,  Nodir  61.5  Deg,  Tan  Alt  40.0  -  54.0  (Km) 
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Scene  5A,  Aug  10.  1996,  LWIR  Nadir  61.7  Deg,  Ton  Alt  57.3  -  71.2  (Km) 
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Figure  3.  False-color  “raw”  images  for  six  MSX  Scenes.  The  scenes  are  obtained  from  the  MSX 
radiometer  “looking”  perpendicular  to  the  satellite  trajectory  at  fixed  nadir  angles  (so-called  “push- 
broom”  scans).  The  “cross-scan”  to  “in-scan”  aspect  ratios  are  much  less  than  1  to  allow  inspection  on 
the  page  (images  appear  compressed  in  the  in-scan  dimension). 
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Figure  4.  Background  radiance  and  relative  standard  deviation  of  radiance  for  six  MSX  scenes.  “Z”  is  distance 
orthogonal  to  the  line-of-sight  and  parallel  to  the  sensor  vertical.  The  average  is  computed  over  the  entire  length  of  the 
scene. 
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Scene  1,  May  25,  1996,  MWIR,  Nadir  1.3  Deg 
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Scene  2,  Nov  13,  1996,  MWIR,  Nodir  55.0  Deg 
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Scene  3,  Aug  17,  1996,  MWIR,  Nadir  57.0  Deg 
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Scene  4.  Aug  19,  1996,  MWIR.  Nodir  61.0  Deg,  Ton  Alt  1.8  -  16.0  (Km) 
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Scene  5.  Auq  10.  1996,  MWIR.  Nadir  61.5  Deq.  Ton  Alt  40.0  -  54.0  (Km) 
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Figure  5.  Segments  of  the  false-color  “detrended”  images  from  Figure  3.  The  images  show  residual 
perturbations  or  radiance  fluctuations  after  removing  the  trends  in  the  in-scan  dimension.  The  images 
have  an  aspect  ratio  =  1 . 
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Figure  6.  In-scan  k a  wave  number  spectra  of  radiance  fluctuations  for  six  MSX  scenes.  The  direction  of  ka  is 
parallel  to  the  sensor  horizontal.  It  is  computed  in  the  direction  of  the  LOS  at  fixed  nadir  angles.  Grey  points  are 
estimates  of  the  averaged  periodograms.  Solid  lines  depict  the  AR  PSD  estimates.  The  straight  line  portions 
show  estimates  of  the  spectral  slopes  in  the  lower  and  higher  wave  number  regions.  The  vertical  bars  depict 
estimates  of  the  standard  error  of  the  averaged  periodogram.  The  instrumental  noise  level  has  been  subtracted 
from  the  AR  PSDs,  but  to  allow  inspection,  has  been  retained  in  the  averaged  periodograms.  The  numbers  in  the 
lower  left  portion  of  each  plot  refer  to  the  image  footprint  (cross-scan  by  in-scan  distances),  TA  =  tangent 
altitude. 
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Figure  7.  Cross-scan  wave  number  spectra  of  radiance  fluctuations  for  six  MSX  scenes.  The  direction  of  £p  is 
parallel  to  the  sensor  vertical  and  is  obtained  orthogonal  to  the  LOS.  Grey  points  are  estimates  of  the  averaged 
periodograms.  Solid  lines  depict  the  AR  PSD  estimates.  The  straight  line  portions  show  estimates  of  the  spectral 
slopes.  The  vertical  bars  depict  estimates  of  the  standard  error  of  the  averaged  periodogram.  The  instrumental 
noise  level  has  been  subtracted  from  the  AR  PSDs,  but  to  allow  inspection,  has  been  retained  in  the  averaged 
periodograms. 
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Spectral  Index  As  Function 
of  Depth  (Numerical  Computation) 
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Simulated  Data  Specifications: 
x,y  Spectral  Indexes  =  -2 
Nx  =  Ny=  1024 
Ax  =  Ay  =  100 

x,y  correlation  lengths  =  1000 
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Figure  8.  Numerical  computation  showing  the  transition 
from  a  k'n  to  a  k(n^1}  spectral  index  as  the  number  of 
correlations  along  the  LOS  increases  in  an  isotropic, 
non-separable  model,  fluctuation  field. 


Figure  9.  Model  calculations  of  the  radiance  variance 
contribution  function  (W(r,0)dr/dz)  vs.  tangent  altitude  for  the 
Bj  band. 


Figure  10.  Model  radiance  spectral  index  as  a  function  of 
local  zenith  angle  with  horizontal  to  vertical  correlation 
length  ratio  =  5. 


Figure  1 1.  Model  radiance  spectral  index  as  a  function  of 
local  zenith  angle  with  horizontal  to  vertical  correlation 
length  ratio  =  20.  Included  on  the  plot  are  solid  square 
symbols  that  show  the  measured  MSX  in-scan  radiance 
spectral  slopes  from  Scenes  1  to  5. 
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